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ABSTRACT 

Fast magnetosonic waves in a magnetically-dominated plasma are investigated. In the 
pulsar wind, these waves may transport a significant fraction of the energy flux. It is 
shown that the nonlinear steepening and subsequent formation of multiple shocks is 
a viable mechanism for the wave dissipation in the pulsar wind. The wave dissipation 
both in the free pulsar wind and beyond the wind termination shock is considered. 
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1 INTRODUCTION 

The mechanisms of the energy transfer from the pulsar to 
the pulsar wind nebula still remain obscure. It is widely ac¬ 
cepted that pulsars emit an electron-positron plasma, which 
form an ultrarelativistic magnetized wind. The rotational 
energy of the neutron star is carried mostly by electromag¬ 
netic fields as Poynting flux (Michel 1982) and should be 
eventually transferred to the radiating particles. The wind 
terminates at a shock front located, in the case of the Crab 
Nebula, some 10^^ cm from the pulsar. The postshock flow 
parameters may be matched with the observed Crab struc¬ 
ture if the energy flux at the shock front is carried mainly by 
the particles (Rees & Gunn 1974; Kennel & Coroniti 1984; 
Emmering & Chevalier 1987). However it is extremely 
difficult for the pulsar wind to reach the kinetic en¬ 
ergy dominated state, at least if the wind is consid¬ 
ered as an ideal MHD flow (Tomimatsu 1994; Bogovalov 
1997; 2001b; Beskin, Kuznetsova & Rafikov 1998; Chiueh, 
Li & Begelman 1998; Lyubarsky & Eichler 2001). The neces¬ 
sary conversion may be provided by dissipation mechanisms. 
These mechanisms operate faster at small scales and first 
affect waves generated in the wind by the rotating, oblique 
pulsar magnetosphere. 

There is a variety of electromagnetic waves in rarefied 
magnetized plasmas (e.g., Akhiezer et al. 1975) however we 
assume that only true MHD waves (those satisfying the con¬ 
dition E = v X B) may be generated by the rotating mag¬ 
netosphere. The reason is that according to the convenient 
view, the plasma density in the pulsar wind is sufficiently 
large such that low-frequency electromagnetic waves are 
heavily damped (e.g., Asseo et al. 1978; Melatos & Melrose 
1996). There are four types of MHD waves but only the 
entropy and the fast magnetosonic (FMS) waves may prop¬ 
agate far beyond the light cylinder where the magnetic field 
is predominantly toroidal. In an oblique magnetosphere, 
the magnetic field at the light cylinder oscillates 
with the pulsar rotation period. These oscillations 


propagate outwards as MHD waves. In the equato¬ 
rial belt of the flow, the magnetic field line at a given 
radius alternates in direction with the frequency of 
rotation, being connected to a different magnetic 
pole every half-period. Such an alternating field is 
transported by the flow in the form of the striped 
wind (Michel 1971, 1982; Bogovalov 1999), which 
may be considered as an entropy wave. At high lat¬ 
itudes, where the magnetic field does not change 
sign, the magnetic oscillations are transferred away 
in the form of FMS waves (generation of FMS waves 
by the rotating, slightly nonaxisymmetric magneto¬ 
sphere was considered by Bogovalov 2001a). 

The entropy wave decays because of the current star¬ 
vation in current sheets separating strips with the oppo¬ 
site magnetic field (Usov 1975; Michel 1982, 1994; Coroniti 
1990). Lyubarsky & Kirk (2001) showed that the flow sig¬ 
nificantly accelerates in the course of reconnection and this 
dilates the timescale over which the wave decays. At typi¬ 
cal conditions, the dissipation radius exceeds the radius of 
the termination shock therefore one should conclude that 
the Poynting flux in the striped wind does not dissipate un¬ 
til the wind enters the termination shock. All the energy 
shonld release within the shock where the flow decelerates. 

In this article, the fate of FMS waves is considered. It 
will be shown that these waves may be described within 
the MHD framework throughout the pulsar wind up to the 
termination shock and ever beyond. The reason is that in a 
magnetically dominated plasma, the FMS waves excite small 
conductivity currents, the oscillations of the magnetic held 
being nearly compensated by the displacement current. In 
MHD regime, FMS waves may decay due to the nonlinear 
steepening and subsequent formation of multiple shocks. It 
will be shown that for typical pulsar parameters, the waves 
may decay only beyond the termination shock. However in 
rapidly spinning pulsars, like the Crab, FMS waves may de- 
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cay before the flow reaches the termination shock provided 
the plasma density in the wind is high enough. 

The article is organized as follows. Properties of FMS 
waves in a magnetically dominated plasma are outlined qual¬ 
itatively in Sect. 2. These estimates are applied to the pulsar 
wind in Sect. 3. The wave decay after the multiple shocks 
formation is considered in Sect. 4. Conclusions are summa¬ 
rized in Sect. 5. In Appendix 1, exact solutions for nonlinear 
FMS waves in the magnetically dominated plasma are pre¬ 
sented. In Appendix 2, the energy and momentum fluxes 
transferred by FMS waves are derived. 


2 FMS WAVES IN A MAGNETICALLY 
DOMINATED PLASMA 


In this section, I consider FMS waves in the wind frame. Let 
us consider waves propagating perpendicularly to the mag¬ 
netic field because far enough from the pulsar, the magnetic 
field is nearly toroidal whereas the waves propagate radi¬ 
ally. Only qualitative estimates are outlined here; the exact 
solutions are presented in Appendix 1. 

When the FMS wave propagates perpendicularly to the 
magnetic field, plasma oscillates along the wave direction 
together with the magnetic field. The flux freezing condition 
ties the plasma density and the magnetic field together: 

— = — = 0 = const. (1) 

n n* 


Here B' and n are the magnetic field and the plasma density 
in the proper plasma frame, the subscript * stands for 
quantities averaged over the wave period. The quan¬ 
tities measured in the proper plasma frame differ 
from those measured in the wind frame because the 
plasma oscillates. In this section, I neglect this dif¬ 
ference assuming the oscillation velocity to be non- 
relativistic. More general consideration is presented 
in Appendix 1. With the aid of Eq.(l), one can express 
the magnetic energy and the pressure via the plasma density 
thus considering the wave as the sound wave in a medium 
with the equation of state 


V =p + 


2 . 

bn 

Stt 


£ = e + 


iZ 2 

bn 

Stt 


where p and e are the plasma pressure and energy den¬ 
sity, correspondingly. The velocity of the wave may be 
now found immediately as the sound velocity in such a 
medium(throughout the paper all velocities are expressed 
in units of the speed of light, c = 1) : 


dV ws^ + 6^n^/47r 
d£ w + tP'n^ 


( 2 ) 


where w = p + e is the enthalpy of the plasma, s = \Jdp/de 
the thermal sound velocity. 

In the pulsar wind, the magnetic energy density signif¬ 
icantly exceeds the plasma energy density. One can conve¬ 
niently express all the values via the parameter 


12 2 
bn 

4:TTW ’ 


( 3 ) 


which is twice the ratio of the magnetic to the plasma energy 
density in the wind frame and the ratio of the Poynting flux 
to the plasma energy flux in the laboratory frame. One can 


see that at ct ^ 1 the FMS velocity is close to the speed of 
light. The corresponding Lorentz factor, 7 m = (1 — 
may be written as 



In the small amplitude, linear wave, the electric field 
may be written as E' = vBi whereas variations in the mag¬ 
netic field are SB' = B'^,{Sn/n.t). Taking into account that 
variations in the density and in the velocity are related by 
the standard expression dn/n* = v/sm, one can see that 

E' = — = SB' (4) 

n* 

to within a factor of 1/a. Because the electric field may not 
exceed the magnetic field, the wave amplitude should be 
limited by 

SB' < B'J2. (5) 


It is shown in the Appendix 1 that when the wave ampli¬ 
tude approaches the limiting value, the oscillation velocity 
becomes ultrarelativistic, the plasma density goes to zero at 
some phase of the wave period and moreover the character¬ 
istic scale for the nonlinear steepening of the wave decreases 
to zero. 

Now let us find what plasma density is necessary to 
allow the MHD solution. The plasma density should satisfy 
the evident condition en > j, where j is the conductivity 
current excited in the wave. Substituting into the Maxwell 
equations 


VxB' = 4^j + —; 

ot 


V X E' = 


as' 


a sinusoidal wave exp(—ioif -I- k • x) with the dispersion law 
ui/k = sm, one can find the conductivity current in this 
wave, 


j = (1 - s\t)iuE'/A- k. 


One can see that in high-cr flows, the conductivity current 
in FMS waves is cr times less than the displacement current. 
Substituting the maximal current density, j = en, and mak¬ 
ing use of Eqs.(2-4), one can see that the cold plasma may 
provide the necessary conductivity current if 

ws > (6) 

where uib = eB'^jm is the Larmour frequency. We consider 
the electron-positron plasma therefore throughout the paper 
m is the electron mass. 

At the conditions (5, 6), FMS waves in a magnetically 
dominated plasma may be considered in MHD approxima¬ 
tion. In the MHD regime, the wave may decay as a result 
of nonlinear steepening and subsequent formation of mul¬ 
tiple shocks. The wave steepening occurs because the wave 
velocity depends on the plasma density and therefore the 
compressive part of the wave moves faster than the expan¬ 
sive one. In the strongly magnetized plasma, even a large 
amplitude wave is only weakly nonlinear because the wave 
velocity varies with the density only by a factor of about 
1/cr. The reason is that, as it was shown above, the con¬ 
ductivity current, which only introduces nonlinearity into 
the MHD equations, is small in the magnetically dominated 
plasma. 
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The difference in the velocities between the compressive 
and expansive parts of the wave may be estimated from 
Eq.(2) as 

^ 1 5n 

OSm ~-• 

cr n* 

The shock forms when the leading edge of the wavefront be¬ 
comes vertical (e.g., Landau & Lifshitz 1959). This occurs 
within a characteristic nonlinear time it takes from the com¬ 
pressive part to shift relative the expansive one by ~ l/u. 
Therefore the shock forms after the wave travels a distance 

*111 ~ ~ u/w. (7) 

This estimate is confirmed by a rigorous derivation in Ap¬ 
pendix 1. After the shock formation, the wave continue to 
distort at the characteristic scale Xni such that eventually 
all the wave energy dissipates in the shock. However both a 
and (u may change because the plasma heats up and acceler¬ 
ates (or decelerates) considerably in the course of the wave 
dissipation. Therefore numerically the dissipation scale may 
differ from the shock formation scale even though the same 
expression (7) is valid for both scales. 


3 FMS WAVES IN THE PULSAR WIND 


Far beyond the light cylinder, the wind may be considered as 
purely radial, whereas the magnetic field as purely toroidal. 
The wind is assumed to be super-FMS, 7 > -s/a. Transform¬ 
ing the electromagnetic fields from the wind frame into the 
pulsar frame, one can see that the condition (5) reduces to 
the condition SB < B* (the factor of two appears because in 
the proper plasma frame = 0 whereas in the pulsar frame 
Et = « B,). So FMS waves may be generated by the 

rotating oblique magnetosphere at high latitudes where the 
magnetic field does not change the sign. In the equatorial 
belt (its width depends on the angle between the magnetic 
and the rotation axes of the pulsar) an entropy wave is gen¬ 
erated in the form of the striped wind. Generally a superpo¬ 
sition of an entropy and FMS waves is generated here (one 
can talk about the superposition because even a large am¬ 
plitude FMS wave is weakly nonlinear in the magnetically 
dominated plasma) however an important point is that an 
entropy wave arises here inevitably because there is no other 
MHD wave that may transfer alternating magnetic field in 
a high-CT plasma. 

Now let us consider validity of MHD approximation 
(Eq.( 6 )) for the FMS wave. The magnetic field in the pul¬ 
sar wind is predominantly toroidal and may be presented 

as (from now and hereafter only averaged quantities 
will be considered therefore the subscript * will be 
omitted) 


( 8 ) 

where Rl = c/Q, = 5 • 10 ®P cm is the light cylinder radius, 
P the pulsar period. The magnetic field at the light cylinder 
may be estimated as 






where = 10®°/iso G-cm® is the magnetic moment of the 
star. The wave frequency in the pulsar frame is the pulsar 


rotation frequency; in the wind frame the frequency is cu = 
17 /( 27 ). Noting that the magnetic field in the wind frame is 
B' — B/ 7 , one can reduce the condition ( 6 ) to 

P < 5-10^'^^cm. (9) 

This radius should be compared with the radius of the stand¬ 
ing shock where the wind terminates. One can place the up¬ 
per limit on this radius by equating the magnetic pressure in 
the wind to the ram pressure of the medium, = pV^, 

where p = Ipo g/cm® is the density of the interstellar gas, 
V — 100V 2 km/s the velocity of pulsar through it. Substi¬ 
tuting Eq.( 8 ), one obtains 

P^hock = 7 X 10^^ cm (10) 

k 2 VP 0 P 

The radius of the terminating shock may be less than (10) 
if the pulsar is surrounded by a dense plerion. Gomparing 
Fq.(lO) with Eq.(9) one can see that FMS waves in the 
pulsar wind may be considered in MHD approximation. 

The wave may decay in multiple shocks that arise as a 
result of nonlinear steepening. Transforming the character¬ 
istic nonlinear scale (7) to the pulsar frame, one gets 

Pni = 47^0-Pz,. ( 11 ) 


The magnetization parameter of the wind depends on the 
plasma density. In the pulsar frame, one can conveniently 
express the density, N = ny, via the dimensionless multi¬ 
plicity factor, K, as 


k,Q,Bl / Rl \ ^ 
27re V P / 


( 12 ) 


Before the multiple shocks arise, the pulsar wind is cold, 
w = mn, and propagates with a constant Lorentz factor 70 . 
The magnetization parameter remains constant, 

B^ COL 

<70 = - = -5 

4TrmNjo 2k7oI7 


where lol = eBLlm is the gyrofrequency at the light cylin¬ 
der. Now one can write the shock formation distance as 

Po = 47 oVoPl = = 2.5 X 10^"^^^ cm. (13) 

KiV kF 

For typical pulsar parameters, 70 ~ 100, k ~ 1 — 10® (Hib- 
schmann & Arons 2001), this radius satisfies the condition 
(9) therefore MHD description of the wave is justified. For 
the typical parameters, this radius exceeds the upper limit 
on the termination shock radius ( 10 ) for all pulsars with the 
exception of the millisecond ones. However there is strong 
evidence to suggest that the multiplicity parameter in the 
Crab pulsar significantly exceeds the typical one and may 
be as high as k ~ 10° (Shklovsky 1970; Rees & Gunn 1974). 
Extra pairs may be produced together with the observed 
gamma radiation in the upper magnetosphere (Cheng, Ho & 
Ruderman 1986; Lyubarskii 1996). At such a plasma density, 
the multiple shocks arise before the flow reaches the termi¬ 
nation shock (observed in the Crab at ~ 10^^ cm from the 
pulsar). So formation of the shocks in the FMS waves may 
occur both in the free pulsar wind within the termination 
shock and in the pulsar wind nebula beyond the termination 
shock. 

If the multiple shocks are not formed until the 
flow reaches the termination shock (i.e., if the ra¬ 
dius (13) exceeds the termination shock radius), the 
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waves enter the nebula through the shock. In case 
the upstream flow is not Poynting dominated, the 
downstream flow is non-relativistic; the nonlinear 
scale (11) is extremely small in this case and there¬ 
fore the waves should decay immediately, in fact al¬ 
ready within the shock. But if the upstream flow is 
still Poynting dominated, the shock is weak and the 
downstream flow is relativistic (Kundt & Krotschek 
1982, Kennel & Coroniti 1984). The waves pass 
freely the shock from upstream because reflection is 
impossible. Parameters of the waves are nearly not 
affected by weak shocks (e.g., Anderson 1963) how¬ 
ever the nonlinear scale decreases because the flow 
decelerates (the statement by Kennel & Coroniti 
that the downstream Lorentz factor remains large, 
~ is based on the assumption that the flow is 
spherically symmetric; one can see that even slight 
deviation of the flow lines from radial may result 
in significant deceleration of the flow). Eventually 
the multiple shocks are formed. Investigation of the 
flow in the nebula and search for the multiple shocks 
formation zone is out of the scope of this article. I 
simply assume that the multiple shocks do form at 
some point and consider evolution of the flow param¬ 
eters in the course of the wave dissipation. It will be 
shown in the next section that in the sub-FMS flow 
(beyond the termination shock) the wave energy dis¬ 
sipates immediately after the multiple shocks arise, 
dissipation being accompanied by the flow deceler¬ 
ation down to 7 ~ 1. Hence although the Poynting 
dominated wind could not be brought to rest in the 
termination shock, the wave dissipation provides an 
effective mechanism of the flow braking. 


4 WAVE DECAY 
4.1 Basic equations 

After the multiple shocks arise, the wave energy dissipates. 
This occurs at the nonlinear scale (11), which now may 
change because the released energy heats the plasma flow. 
One can consider evolution of the flow parameters in the 
conrse of the wave decay applying the conservation laws to 
the system. Let us separate the average fluxes of conserv¬ 
ing values into the wave and the flow parts. Such a separa¬ 
tion may be performed in general terms for small amplitude 
waves (see Appendix 2). The average plasma density, n, and 
velocity, v, (and the corresponding Lorentz factor 7 ) are de¬ 
fined such that the wave does not contribnte to the matter 
flux (in Appendix 2 the averaged quantities are marked by 
tilde; here we work only with averaged quantities and tilde is 
omitted). Therefore the continuity equation may be written 
as (the flow is considered as locally spherically symmetric) 

njvR^ = myovoRo, (14) 

where the subscript ” 0 ” is referred to the quantities at the 
shock formation radius (13). The energy equation may be 
written as 


where Q is the energy transferred from the wave to unit 
plasma volume per unit time. The magnetic field strength 
may be expressed via the density making use of the frozen 
flux condition (see Eq.(A1.16)) 


B ^ Bq 
R ny Ronoyo ’ 


(16) 


The entropy equation describes the entropy growth as a re¬ 
sult of the wave decay: 


dr 


d 1 
- 

dr n 


+ 3 


.dT 
dr ’ 


(17) 


where dr = dR/y is the proper time. We assume that just 
after the shock formation, the plasma is heated to relativistic 
temperatures, w = dnT. 

Because the wave propagates in the proper plasma 
frame with the velocity sm = 1 — 0 (l/cr), its energy and 
momentum are equal, in this frame, to within 1/a (see Ap¬ 
pendix 2). Therefore if the energy de per unit mass is ab¬ 
sorbed in the proper plasma frame, then the momentum 
de per unit mass is also transferred from the wave to the 
plasma. Then in the pulsar frame, the absorbed energy is 
dE = 2yde. Taking into account that dVdt is invariant, one 
can write the energy release rate as 


Q 


dE 

dVdt 


2yn 


de 

dr 


Substituting de = TdS, one writes the entropy equation in 
the form 



„ dT 
^^IR 



(18) 


Let us now turn to the energy equation (15). Making 
nse of Eqs.(14, 16), one can see that the Poynting flux (the 
second term in the brackets in the left-hand side) is nearly 
independent of R in the ultrarelativistic case. Expanding 
this term in I /7 to the next order and retaining in the first 
term only leading terms in I /7 because this term is already 
small as 1/a, one can reduce Eq.(15) to 


d 

'dR 



BlRl „2 

STry^ 


Qi?^ 


(19) 


Substituting n from Eq.(14) into Eqs.(18, 19), one gets 
finally 

+ QR" . ( 20 ) 

dR\Toyo 372 y mylnoRl' ^ ’ 

7 d{T/To) T d{y/yo) 2Ty QR^ 

70 dR To dR ToyoR 2Tono7o.Ro' 

where the factor 


SBq _ 3fro 

2,2-KnoToyl 2yl' 


( 22 ) 


is less than unity in super-FMS flows and exceeds unity in 
the opposite case. It was found above that at different wind 
parameters, the multiple shocks may arise in FMS waves 
both before the wind reaches the termination shock and be¬ 
yond the shock in the pulsar wind nebula. Therefore one 
should consider the wave decay both in the super-FMS (free 
wind) and in the sub-FMS (beyond the termination shock) 
flows. The flow behavior in these cases is just opposite there¬ 
fore let us consider them separately. 
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4.2 Wave decay in the super-FMS flow 


where 


The flow in the wind is super-FMS, in this case a ^ 1 and 
the plasma is cold at Ro so one can take To ~ m. Eliminating 
Q from Eqs.(20, 21) and neglecting the term with a, one gets 

7 d{T/To) T d{y/yo) 2Tj 

70 dR To dR To^oR ’ 

which yields 

rji r?2 rp p2 

W . (23) 

7 70 

So in the super-EMS flow heating of the medium is accom¬ 
panied by acceleration. 

The energy release is caused by the nonlinear wave dis¬ 
tortion. Therefore the fraction of the wave energy dissipated 
at the distance Ai? may be roughly estimated as A7?/i?ni. 
One should substitute in the expression for the nonlinear 
scale ( 11 ) parameters varying in the course of the energy dis¬ 
sipation. Variation of a may be found making use of Eqs.(14, 
16, 23): 

(iorV 

16'KnT'y \^7i?o J 


Jro 4noTo7g 

is the ratio of the transferred energy to the initial plasma 
energy. When the wave decays, x becomes very large and 
reaches the initial ratio of the wave energy to the plasma 
energy; for large amplitude waves this ratio is about no- 
Integrating Eq.(24), one yields 


^7 

To7o 


= l + -+x- 


a7o 
87 ^ ' 


Substituting T from this relation into Eq.(25), one gets the 
equation 


2 

7 


1 +-+*- 



d') 

dx 


1 , 


which is linear with respect to *( 7 ). The solution is 


7g 


3 -I- a -I- 3* ± ■y/9(l — a)2 -|- 6x(3 -I- a) -I- 
2(3 - a) 


(26) 


Now the dissipated fraction of the wave energy may be es¬ 
timated as ~ AR/{a'y^R l) = RoAR/R?. One can see that 
the significant part of the wave energy dissipates at the scale 
AR ~ Ro - The dissipation scale is roughly equal to the shock 
formation scale because acceleration of the flow is compen¬ 
sated by the plasma heating and corresponding decreasing 
of G. 

Taking into account that the wave energy initially ex¬ 
ceeds the plasma energy about go times (for large ampli¬ 
tude waves), one can see from Eq.(20) (at a ^ 1) that 
T 7 /(To 7 o) ~ Go after the wave dissipates. This implies, 
with the aid of Eq.(23), that T ~ m^/Go{Ro/R), 7 ~ 
y/Go'yo{R/Ro)- So the wave dissipation at 7? ~ i?o heats 
the plasma till the temperature ~ m^/Go and accelerates 
it to the Lorentz factor ~ v^^ 7 o- Then the heated plasma 
cools and accelerates in the Bernoulli regime, T oc 1/R, 
7 oc 7?. At the distance Ri ~ y/aoRo, all the released energy 
transforms to the kinetic energy and the plasma reaches the 
Lorentz factor 


71 


<7070 


2kQ 


= 1.3 X 10 


7 U30 

kP^' 


4.3 Wave decay in the sub-FMS flow 


If the multiple shocks formation radius (13) exceeds the ter¬ 
mination shock radius the waves pass the termination shock 
and enter the sub-FMS flow. In the sub-FMS case, the flow 
decelerates in the course of energy release (see below) and 
the nonlinear scale sharply decreases just after the multiple 
shocks formation. Therefore the wave decay scale turns out 
to be small as compared with the radius and one can con¬ 
sider the flow in the plane geometry. In this case Eqs.(20, 
21 ) reduce to 


— ( ^ 1. 
dx \Toyo J 


(24) 


3 7 d{T/To) T < 7 ( 7 / 70 ) 

70 dx To dx 


(25) 


where the sign of the square root should be positive for a < 1 
and, correspondingly, negative for a > 1 to satisfy the initial 
conditions. One can see that the super-FMS flow (a < 1) 
accelerates whereas the sub-FMS flow (a > 1) decelerates 
in the course of energy release. For a > 1, x ^ 1 one finds, 
with the aid of Eq.(22), 



It follows from Eqs.(14, 16) that in the narrow decay zone, 
R — Ro ^ Ro, the magnetic field remains constant while 
the flow remains relativistic, u « 1. The absorbed energy is 
transformed into the internal plasma energy and the mag¬ 
netization parameter decreases: 

= B^ _ ^ 

4Trwj^ X 

The wave energy dissipates completely when a; ~ < 7 o S> 1. 
The nonlinear scale (11) decreases already when x exceeds 
unity i.e. when the fraction of the dissipated energy is small 
(~ l/<7o) therefore the total dissipation scale turns out to 
be small, ~ Ro/go- 

When the wave decays completely, x ^ go and one gets 
7 ~ 1, (7 ~ 1. However this solution becomes invalid when, 
in the proper plasma frame, the particle Larmour radius, 
Tg = r/(e73'), exceeds the wavelength, A = 2n'yRL- Making 
use of Eqs.( 8 , 12, 22, 27 ), one can write the condition Vg ~ 
2'k^Rl as 

X /‘i2'KKRL\'^^^ 

Go V Ro J 

The right-hand side of this expression is small therefore only 
a small fraction of the energy is transferred to the plasma 
in MHD regime; most of the energy is dissipated when the 
Larmour radius exceeds the wavelength. This may result 
in formation of a high energy tail in the particle energy 
distribution. 
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5 CONCLUSIONS 

In the pulsar wind, a significant fraction of the Poynting flux 
may be transported by FMS waves. The plasma in the pul¬ 
sar wind is magnetically dominated in the sense that in the 
proper plasma frame, the magnetic field energy density ex¬ 
ceeds the plasma energy density. This condition is equivalent 
to the condition that in the pulsar frame, the Poynting flux 
dominates the plasma energy flux. In such a plasma, FMS 
waves excite small conductivity current, oscillations of the 
magnetic field being nearly compensated by the displace¬ 
ment current. Thus a very small plasma density is sufficient 
to keep this current and MHD description of these waves 
is justified throughout the pulsar wind till the termination 
shock and beyond. 

These waves may decay in multiple shocks that arise 
through nonlinear steepening of the waves. Because the con¬ 
ductivity current is small at the pulsar wind conditions, even 
a large amplitude FMS wave is nearly linear and therefore 
shocks are formed at large distances from the pulsar. De¬ 
pending on the plasma density in the wind, this may 
occur either before the flow reaches the termination 
shock or beyond the shock. In any case the wave 
energy adequately dissipates in the multiple shocks. 
Energy release in the super-FMS flow is accompa¬ 
nied by the plasma heating and acceleration. The 
thermal pressure also does work on the flow and 
the plasma accelerates further on. Therefore if the 
waves dissipate in the free, super-FMS pulsar wind, 
all the wave energy eventually transforms into the 
kinetic energy of the flow. If the waves do not dis¬ 
sipate in the wind, the flow remains Poynting dom¬ 
inated. In this case the termination shock should 
be weak and the downstream flow remains ultra- 
relativistic. Such a flow does not match the slow 
expansion of the nebula; this was considered as an 
evidence for conversion of a significant fraction of 
the Poynting flux into the particle energy flux before 
the pulsar wind reaches the termination shock (Rees 

6 Gunn 1974; Kennel & Coroniti 1984). However it 
was shown above that the wave dissipation in the 
sub-FMS downstream flow is accompanied by the 
flow deceleration to subrelativistic velocities. Hence 
even if the flow remains Pointing dominated at the 
termination shock, the wave dissipation downstream 
the shock may provide the necessary deceleration. 

Only a fraction of the total Poynting flux is transferred 
by the FMS waves therefore the proposed mechanism does 
not solve the u-problem but should be considered as an el¬ 
ement of the futnre complete theory. The mechanism con¬ 
sidered may also play an essential role in gamma-ray burst 
models involving Poynting dominated outflows from com¬ 
pact objects (Usov 1992, 1994; Meszaros & Rees 1997; Kluz- 
niak & Ruderman 1997; Blackman & Yi 1998; Spruit 1999; 
Lyutikov & Blackman 2001; Drenkhahn 2002; Drenkhahn & 
Sprnit 2002). 
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APPENDIX 1. NONLINEAR FMS WAVES IN A 
MAGNETICALLY DOMINATED PLASMA 


Relativistic nonlinear MHD waves are considered in the gen¬ 
eral case by Akhiezer et al (1975). Here we consider only the 
waves propagating perpendicularly to the ambient magnetic 
field in the case a 1. 

Let us first consider the waves in the plane geometry. It 
is followed from the continuity equation 


d d ^ 

-^'rn + = 0 


(Al.l) 


and the frozen-in condition 


dB d ^ ^ 
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that the magnetic field may be presented in the form 

B = bn-y, b = —= const, (A1.2) 

n*7* 


where the subscript * stands for quantities at some 
fiducial point; one can conveniently choose this point 
in such a way that B* is equal to B averaged over 
the wave period. The dynamic equation may be presented 
in the form of the energy equation 


dTw dToi 
dt dx 


(A1.3) 


where the components of the energy-momentum tensor are 


rp 2 . B 2 

ioi = wy V + -—y V, 
Atv 


(AlA) 


Too = wy^ —p + ^ ^ ^ B^y^. (yll.5) 

OTT 

Note that this set of equations is reduced, with the aid of 
Eq.(A1.2), to the standard hydrodynamics equations with 


E = e + 


Stt ’ 


V ^ p + 


_ 

Stt 


Kn^ + 


Fn\ 
Stt ’ 


d 2 

\p = S + 'P^e+p+—. (A1.6) 

47r 

The nonlinear simple wave may be found from the con¬ 
dition that all dependent variables are functions of one of 
them, e.g., n (Landau & Lifshitz 1959). This means that 
Eq.(A1.3) should be equivalent to Eq.(Al.l), i.e. 


dToi _ dinyv) 
dToo d{ny) 


The last equation reduces to 


2 dv sm 

T' = -’ 

dn n 


(A1.7) 


(A1.8) 


where the FMS velocity sm is given by Eq.(2). For the small 
wave amplitude, Sn = n — n* <C n,, one obtains the linear 
FMS wave propagating with the phase velocity sm = const. 
For a strong wave, Sm depends on the local density therefore 
the wave becomes nonlinear. However in the strongly mag¬ 
netized case, 0^1, Sm goes to unity and therefore even a 
strong wave is nearly linear. 

Substituting sm = 1 into Eq.(A1.8), one obtains 


n_ 1 V 1 — V* 
n* y I — V I + V* 


(A1.9) 


In the frame moving, in average, with the plasma, u* = 0, 
one gets for the small amplitude wave 


V = dn/n*. 

In the case 7 1 Eq.(A1.9) reduces to 

7 = 7*n/n,. 


The waveform moves along the characteristic of 
Eq.(Al.l): 

dx _ d{nyv) _ n + sm _ 2'Kwnl{l — 1 — u* 

dt d{ny) 1 -|- vsm b'^n* 1 -I- w* 

(Al.lO) 


The last equality in this expression is obtained for the high¬ 
er case. Let us consider the cold plasma, w = mn. Then 
Eq.(Al.lO) yields 



1 1 — u* 
2 (t* 1 -|- u* 



t + f{n), 


(Al.ll) 


where cr, = 47rm/(fe^n*). The function /(n) is determined 
by the initial waveform. For the initially sinusoidal wave. 


B = B* (1 -I- acoscjx) 
me gets from Eqs(A1.2, A1.9): 


f(n) = — 


1-1-u* \( n\'^ 



(A1.12) 


Note that the dimensionless amplitude, a, does not exceed 
the value Umax = (1 -I- u*)/2; n ^ 0 at some phase of the 
wave period as a ^ cimax ■ 

According to Eq.(Al.lO), any portion of the wave ad¬ 
vances with a constant velocity depending only on the 
plasma density. The compressive parts overtake the expan¬ 
sive ones and eventually a shock arises where the waveform 
becomes vertical. 


dx d?x 
dn dn'^ 


This occurs at a time 


to — 


/To 20-, 1 -h 

Y 3 Qlo 1 — Vt 


^4 — y/l + lb{a/amaxY'^ 
y/l + - 1 

(A1.12) 


When the wave amplitude approaches the maximal one, the 
shock formation time goes to zero; at small amplitudes, to 
grows asl/a (Fig. 1). At intermediate amplitudes, the shock 
formation time is to ~ cr*/a; in the frame moving with the 
plasma, u* = 0; if the plasma moves with respect to the 
observer with an ultrarelativistic velocity, 7* ;?> 1, the shock 
formation time is to ~ atylfio. 

Now let us consider the radial flow. In this case the 
governing equations are written as 




(A1.13) 


^Too + Toi = 0 . (A1.15) 

It follows from Eqs.(A1.13, A1.14) that one can present the 
magnetic field as 


B = hRny, 


(A1.16) 


where b = const. If the flow is cold, w = nm, Eqs.(A1.13, 
A1.15, A1.16) reduce to Eqs. (Al.l, A1.2, A1.3) by substi¬ 
tution n — nB?. Therefore the above estimate of the shock 
formation time remains valid also for the radial flow. 


APPENDIX 2. ENERGY AND MOMENTUM OF 
THE WAVE 

Let us show that for small amplitude waves, one can separate 
the energy and momentum of the wave and the flow without 
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ala 


max 


order terms expanding Eqs.(A2.1, A2.2) in 5v and 5n to the 
second order. The result is 

< Sn >= < {Sv)^ >; 

< 6v >= —— < > . 

sm 

Now expanding Eqs.(A2.3, A2.4) and making use of the ther¬ 
modynamical expression dS/dn = Wln, (becanse the wave 
is isentropic, all thermodynamical values may be considered 
as functions of n), one yields 

< Too >= -h W* < >; 


< Toi >= smW. < {5vf >, 

where = £{n,) etc. So the energy density is separated 
into a part which depends only on the average medium den¬ 
sity and the wave part which is proportional to the wave 
amplitude squared. The momentum density of the medium 
is zero in the proper plasma frame therefore only the wave 
momentum contributes to < Toi >• The energy and mo¬ 
mentum of the wave coincide, in the proper plasma frame, 
to within 1 — sm Ijcr. Making the Lorentz transform, one 
can get components of the energy-momentum tensor in an 
arbitrary frame of reference. 


Figure 1. Dependence of the shock formation time, to, on the 
wave amplitude, a. 


considering specific waveforms. Let us define the average 
plasma density, n*, and velocity, «*, by the relations 

< nwy >= n*w*7*; (^2.1) 

< ny >= n*7*, (A2.2) 

where the angular brackets denote averaging over the wave 
period. With such a definition, there is no mass flux associ¬ 
ated with the wave. 

The total energy and momentum densities are (see 
Eq.(A1.4-A1.6)) 

Too = - V; (A2.3) 

Toi = Wvy^. (A2.4) 

Averaging this values, one can present them as a superpo¬ 
sition of a flow part dependent only on average plasma pa¬ 
rameters and a wave part dependent on the wave amplitude. 
Calculations are simplified in the frame where plasma is at 
rest in average, v* = 0. 

In a small amplitude wave, the energy and momentum 
are of the second order in the wave amplitude therefore lin¬ 
ear relations, like (see Eq.(A1.8)) 

Sn 5v 
n* Sm ’ 

may be used only in the second order terms. The average 
of the first order terms may be expressed via the second 
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